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Abstract 

We solve the spectrum of quantum spin chains based on representations of the Temperley- 
Lieb algebra associated with the quantum groups U q (X n ) for X n = A\, B n , C n and D n . The 
tool is a modified version of the coordinate Bethe Ansatz through a suitable choice of the Bethe 
states which give to all models the same status relative to their diagonalization. All these 
models have equivalent spectra up to degeneracies and the spectra of the lower dimensional 
representations are contained in the higher-dimensional ones. Periodic boundary conditions, 
free boundary conditions and closed non-local boundary conditions are considered. Periodic 
boundary conditions, unlike free boundary conditions, break quantum group invariance. For 
closed non-local cases the models are quantum group invariant as well as periodic in a certain 
sense. 



1 Introduction 



Temperley-Lieb (TL) algebra Jl[] has been widely used in the construction of solutions of the Yang-Baxter 
equation [^], |3|], which is a sufficient condition for integrability of lattice models through the quantum 
inverse scattering method. On the other hand, quantum group[f§ in some cases is a symmetry of the 
integrable model [Efl. Thus, quantum group together with Temperley-Lieb algebra provide a powerful 
algebraic framework to build and study various kinds of lattice models in two-dimensional statistical 
mechanics. One particular way of building models which are quantum group invariant uses the TL 
algebra. It is a algebra generated by the Hamiltonian density U k , k = 1,2, ...,N — 1 subject to the 
following constraints 

U 2 k = (Q + Q- l )U k , U k U k±l U k = U k 

(1.1) 

UkUj^UjUk \k-j\>l. 

with Q E C a given number. This algebra appears in a large class of quantum lattice models and leads, 
at level of free and ground state energies, to some equivalence among the models 0. 
Taking into account usual boundary conditions, the TL Hamiltonians take the form 

N-l 

n= J2 u kk+i + kt (i.2) 
fe=i 

where U kk+ i = U k operates in a direct product of complex spaces at positions k and k + 1. In general, 
they are not invariant with respect to quantum groups since the boundary terms, bt, break translational 
invariance, reflecting the non-cocommutativity of the co-product. Indeed, we know from |4], [7] that very 
special boundary terms must be considered when we seek these quantum group invariant spin chains. In 
particular, one possibility is to consider free boundary conditions, i.e., bt = . For The XXZ-Hamiltonian 
with free boundary conditions one has to apply the Bethe ansatz || techniques introduced first by Alcaraz 
et al and after by Sklyanin [|] using Cherednik's reflection matrices |TD|, [11]. By this method the 



XXZ-Heisenberg model |[12|| , the spl q (2,l) invariant supersymmetric t-J model []T3|, , the U q [sl(n)] 
invariant generalization of the XXZ-chain []T5[] and the SU q (n\m) spin chains []7|, |HJ have been solved for 
free boundary conditions. 



Recently, by means of a generalized algebraic nested Bethe ansatz, Karowski and Zapletal JT7[] presented 



a class of quantum group invariant n-state vertex models with periodic boundary conditions. Also an 
extension of this method to the case of graded vertex models was analyzed in Jl8|], where a spl q (2\l) 
invariant susy t-J model with closed boundary conditions was presented. 

In fact, these kind of models were first discussed by Martin from the representations of the Hecke 
algebra. In this case, the boundary term is a non-local operator, bt = U (see Section 4). The study of 
closed quantum group invariant closed spin chains in the framework of the coordinate Bethe ansatz was 
presented by Grosse et al [T^DJ for the fundamental representation of SU q (2) and generalized by [|21[| . In 



this context it would be interesting to discuss other quantum group invariant closed spin chains. 

In the present paper we find the spectrum of Hamiltonians based on representations of the TL algebra 
associated with quantum groups. The Bethe ansatz equations for different types of boundary conditions 
(periodic, closed and free) are obtained through a modified version of the coordinate Bethe ansatz. 
Therefore we generalize the results of ref.^lJ.Our method was recently applied to solve graded T-L 



Hamiltonian [^2] and anisotropic correlated electron model associated with this algebra 

The paper is organized as follows. In Section 2, we describe the representations of the TL algebra, 
constructed as projectors on total spin zero of two neighboring spins. In Section 3, we introduce the 
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modified coordinate Bethe Ansatz for the TL Hamiltonians with periodic boundary conditions. In Section 
4 their Bethe Ansatz solution is presented with non-local boundary conditions. In section 5, contains the 
solution for free boundary conditions. Finally the conclusions are reserved for section 6. 



2 Representations of the TL algebra 

Representations of the TL algebra, commuting with quantum groups, can be constructed in the following 
way ||24{|. Suppose U q (X n ) is the universal enveloping algebra of a finite dimensional Lie algebra X n , 



equipped with the coproduct A : U q — > U q ® U q p5[. If now tt : U q — > End V\ is a finite dimensional 



irreducible representation with highest weight A and we assume that the decomposition V\ <8> V\ is 
multiplicity free and includes one trivial representation on Vq, then the projector Vq from V\ (g> Va onto 
Vq is a representation of the TL algebra. The deformation parameter q, which plays the role of a coupling 
constant in the Hamiltonian, is related to Q as: 

Q + Q- 1 = Tr VA (<T ap ), (2-1) 

where p is half the sum of the positive roots. 

We will consider the following specific cases, (V\,U q (X n )) = (V 2s A 1 ,U q (Ai)) for sp\n s, (V/ Vl ,U q (B n ),(V^ 1 ,U q (C 
and (V\ 1 ,U q (D n ). Namely we treat the g-deformations of the spin-s representation of sl(2) and the vector 
representations of so(2n+ l),sp(2n) and so(2n). V\ denotes the U q (X n ) module with highest weight A 
and Ai is a highest weight of X n . 

Since we are not going to use any group-theoretical machinery, we will just lift the relevant formulas 
off Batchelor and Kuniba in order to display explicitly the Hamiltonians to be diagonalized. 

We introduce orthonormal vectors e M , < e M , e v >= 5^, to express p and the set A of weights appearing 
in the representation n of U q (X n ) as follows: 

A\ : Spin-s representation of the sl{2) algebra 

A = {sei2, (s - l)ei2, . . . , -se 12 }, p = -e 12 , (e 12 = e 1 -e 2 ), 
J={s, 8 -l,...,-s}, e(/i) = (-lf, Q + Q- 1 = [2s + \]. (2.2) 

B n : Vector representation of the so(2n + 1) algebra (n > 2) 

n i 

A = {0, ±ei, . . . , ±e n }, p=Y J (n-(a- -))e Q , J = {0, ±1, ... , ±n}, 

«M = (-!)', Q + Q- ^ 2 "- n y . (2.3) 
C n : Vector representation of the sp(2n) algebra (n > 1) 

n 

A = {±ei, . . . , ±e n }, p = ^2(n - a + l)e a , J = {±1, . . . , ±n}, 

a=l 

c(a*) = signM, Q + Q- 1 = [n ] [2n + 2] ■ (2.4) 

[n + 1J 
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D n : Vector representations of the so(2n) algebra (n > 3) 



n— 1 



A= {±e 1} ...,±e n }, p=J2( n ~ a ) e a, J = {±l,...,±n}, 

0=1 

e 00 = l, Q + Q-=!?p«. (2.5) 

[n - 1J 

For ll £ J the symbol fx is defined as fx — fx + 1/2 for Ai with s semi-integer and fx = fx for A\ with 
s integer. For _B n , C n and D n , fx = with the exception of = 1 for B n . The g-number notation is 
[x] = (q x — q~ x ) /(q — q~ r )- 

For X n = B n ,C n and D n , we extend the suffix of to — n < ll < n by setting e_ M = — e M (hence 
e = 0). Using the index set J above, we can write A = {/x(ei — e 2 )} for A 1 and A = {e^|/i G J)} for 
5 n ,C„ and £) n . 

Denoting by e EikIVa the matrix unit, having all elements zero, except at row ll and column z/, 
the projector can be written as 

V ° = TTTTvT E e^e^q-^+^^E^ ® (2.6) 
v + v Mi i, 6 j 

In the following we will refer to all models generically as higher spin models for simplicity, even when 
not talking about A\. The dimensions of B n are related with the dimensions of A 1 for s integer and the 
dimensions of C n and D n are related with the dimensions of A 1 for s semi-integer. 

If we consider then a one-dimensional chain of length N with a "spin" at each site, the "spin variables 
range over the set of weight vectors A = {v^\ll e J} and our Hilbert space is an iV-fold tensor product 
Va <S> ■ ■ ■ <8> V\. For A\, these are the g-analogs of the usual spin states. The periodic Hamiltonians 
associated with the TL representations are given by the following sum over N sites 

N 

H = J2U k . (2.7) 
k=i 

Here bt = Un n+i — Uni and the Hamiltonian densities acting on two neighboring sites are then given 
by: 

(v,\\U\h,k) = e(/i)e(z/)g- <e - +e -^^ +K , ^ + A,o 

fX,V,K,X e J (2.8) 

Using representations of the TL algebra, one can also build solvable vertex models whose Hamiltonian 
limit leads to the previously mentioned quantum spin chains. To do so, we introduce an operator R(u) E 
End (Va <g> V A ) by 

R k (u) = I <g> • • • <g> I <g> R(u) <g> I <g> ■ ■ ■ ® / 

fc,fe+i 

*(«) = E Rfc^E™ ® 
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where u is the spectral parameter and the anisotropy parameter r\ is chosen so that 

2 cosh rj = Q + Q' 1 (2.10) 
The R matrix commutes with the quantum group action and the Yang-Baxter equation 

R k (u)R k+1 (u + v)R k {v) = R k+l (v)R k (u + v)R k+1 (u) (2.11) 

is valid owing to the TL relations 

Due to the TL algebra these models are equivalent to the 6-vertex and the " p-state " self-dual Potts 
models ( Jp = Yim q ^i(Q + Q^ 1 )) through the argument in 0. In fact, the cases U q {Ai) with s = 1/2 
and U q (Ci) yield the 6-vertex model itself. When q = 1, the vertex models here reduce to those discussed 
in [^6], [2~8|] , where the number of states is equal to v /p. 

The case A\ have been studied by several authors. When q = 1 and s = 1/2 the operator U k 
is essentially the Heisenberg interaction term cr^cr k+1 + <^ k a k+1 + cr k cr k+1 . Besides the explicit matrix 
elements , available from (|2~%|), the local Hamiltonian U k is in principle also expressible in terms of the 
usual representation matrices for SU(2) generators. The resulting g-deformed Hamiltonian has been 
written down for s = 1 in 0. However for s > 1, writing the Hamiltonian in terms of the usual SU(2) 
operators becomes very cumbersome. These Hamiltonians have alternative expressions in terms of Casimir 
operators 0. 



The limit q — > 1 has discussed for general s p9[ [30|1 . As for related spin-1 biquadratic model [^9], ^lj] 



they are massive for s > 1 and of relevance to the dimerization transition on SU(n) antiferromagnetic 
chains The case s = 1/2 has also been investigated in some detail in [ffl]. 



Having now built common ground for all models, whose salient feature is they being spin zero projectors, 



we may now follow the steps of reference to find their spectra 



3 Bethe Ansatz: Periodic boundary conditions 

All the above Hamiltonians are U(l) invariant and we can classify their spectra according to sectors. For 
A 1 the commuting operator is the total spin S z = J2k=i $k where S k =diag(s, s — 1, — s + 1, — s) k and 
we set the conserved quantum number r = sN — S z . We extend these quantum numbers for the other 
algebras as r = Nlo — S z where S z = J2k S| with S z k =diag(tu, uo — 1, —lu + 1, —uo) k and u = max{ J}. 
Eigenvalues of the operator r can be used to collect the eigenstates of H in sectors, \l/ r . Therefore, there 
exists a reference state satisfying H^o = E ^ , with E = 0. We take \I>o to be = \lo uo lo ■ ■ ■ uo) . 
It is the only eigenstate in the sector r = and all other energies will be measured relative to this state. 

In every sector r there are eigenstates degenerate with \I>o- They contain a set of impurities. We call 
impurity any state obtained by lowering some of the \u, k)'s, such that the sum of any two neighboring 
spins is non-zero. Since TC is a projector on spin zero, all these states are annihilated by H. In particular, 
they do not move under the action of H, which is the reason for their name. 

Nothing interesting happens in sectors with r < 2uo. Sector r = 2uo, we encounter the situation where 
the states \a,k) and \—a,k±l), a 6 J, occur in neighboring pairs. They move under the action of 
H, i.e., the sector r = 2uo contains one free pseudoparticle. In general, for a sector r we may have p 
pseudoparticles and iV w _i, A^_ 2 , N^+i impurities of the type lo — 1, u — 2, — uo + 1, respectively, 
such that 

2cj-1 

r = 2uop+ <*N u -a. (3.1) 

a=l 
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The main result of this section is to show that H can be diagonalized in a convenient basis, constructed 
from products of single pseudoparticle wavefunctions. The energy eigenvalues will be parametrized as a 
sum of single pseudoparticle contributions. 

The first nontrivial sector is r = 2uj and the correspondent eigenspace is spanned by the states 

\k(—a, a) >= \uj to ■ ■ ■ to -a a uj ■ ■ ■ uj > , where k = 1, 2, N — 1 and a E J. We seek eigenstates of 

k 

Ti which are linear combinations of these vectors. It is very convenient to consider the linear combination 

UJ 

\Q(k)) = £ e(u)e(a) q ~<^P> \k(-a, a)) , (3.2) 

a=—ui 

which is an eigenstate of 

u k \n(k)) = (Q + Q- 1 )\n(k)). (3.3) 

Moreover, the action of Uk±i on \Q(k)) is very simple 

u k - x \n(k)) = e \n(k - 1)) u k+ i \n(k)) = e \n(k + 1)> 

(3.4) 

U k \n(m))=0 k^{m±l,m} 

where e = 1 for B n , D n and Ai ( s integer) and e = — 1 for C n and Ai ( s semi-integer). 

It should be emphasized that although the linear combination ( |3.2| ) is different for each model, the 
action of U k is always given by ( p.3| ) and (|3.4|). Therefore, all Hamiltonians can be treated in a similar way 
and it affords a considerable simplification in their diagonalizations when we compare with the calculus 
used in the usual spin basis [B3] . 



3.1 One-pseudoparticle 

We will now start to diagonalize H in every sector. Let us consider one free pseudoparticle as a highest 
weight state which lies in the sector r = 2u 

^ = E^)l fi (*0>- (3-5) 

k 

Using the eigenvalue equation Ti. ty 2 u> = E 2ul ^2u, one can derive a complete set of equations for the 
wavefunctions A{k). 

When the bulk of H acts on \£l(k)) it sees the reference configuration, except in the vicinity of k 
where we use ( |3~3"| ) and (|3.4| ) to get 

H \Vt{k)) = (Q + Q- 1 ) \tt(k)) + e p(k - 1)) + e \tt(k + 1)) 

2<k<N-2 (3.6) 

Substituting (|3~6| ) in the eigenvalue equation, we get 

(E 2ul -Q- Q- l )A(k) = eA(k - 1) + eA(k + 1) 

2<k<N-2 (3.7) 

Here we will treat periodic boundary conditions . They demand Un,n+i = Un,i, implying A{k + N) = 
A(k). This permits us to complete the set of equations ( |3.7| ) for A(k) by including the equations for 
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k = 1 and k — N — 1. Now we parametrize A(k) by plane wave A{k) = A^ k to get the energy of one 
free pseudoparticle as: 

E 2w = Q + Q~ 1 + e(£ + C 1 ) 

e = 1 (3.8) 

Here £ = e ld , 9 being the momenta determined from the periodic boundary to be 9 = 2ttI/N, with / an 
integer. 

3.2 One-pseudoparticle and impurities 

Let us consider the state with one pseudoparticle and one impurity of type (a> — 1), which lies in the sector 
r = 2u> + 1. We seek eigenstates in the form 

^+i(6,6)= E {A 1 (k 1 ,k 2 )\n 1 (k l ,k 2 )) + A 2 (k 1 ,k 2 )\n 2 (k 1 ,k 2 )}} (3.9) 

fci<A;2 

We try to build these eigenstates out of translational invariant products of one pseudoparticle excitation 
with parameter £ 2 and one impurity with parameter ^: 

*Wi(£i,6) = |(w - l)(£i)) x * 2 ,(6) + * 2 .(6) x |(w - 1)^)) (3.10) 
Using one-pseudopaticle eigenstate solution ( |3.5| ) and comparing this with ( |3.9| ) we get 

k 2 )) = £ e(w)e(a) g-< e " +e «^> \k x {u - 1), fc 2 (-a, a)) 
|fi 2 (*i, fc 2 )) = £ e(w)e(a) g-< e "+ e -"> ^(-a, a), fc 2 (cu - 1)) 

(3.11) 

and 

A 1 (k ll k 2 ) = A^ 2 , A 2 (h, k 2 ) = A^tf 2 - (3.12) 

Periodic boundary conditions A\(k 2 , N+ki) = A 2 (ki, k 2 ) and Ai(N + k\, N + k 2 ) = Ai(ki, k 2 ), i — 1, 2 
imply that 

4^ = 4, , ^ = (66)^ = 1 (3.13) 

When H now acts on ^ 2uJ+ i, we will get a set of coupled equations for Ai(k\, k 2 ), i = 1, 2. We split 
the equations into far equations, when the pseudoparticle do not meet the impurity and near equations, 
containing terms when they are neighbors. 

Since the impurity is annihilated by H, the action of H on ( |3~9"1) in the case far (i.e., (k 2 — k{) > 3), 
can be write down directly from ( |3.7|) : 

(e 2uj+1 - Q- Q- 1 ) A^h, k 2 ) = eA 1 (k 1 , k 2 - 1) + eAxih, k 2 + 1) (3.14) 
(E^+i - Q- Q- 1 ) A 2 (k h k 2 ) = eA 2 (h - 1, k 2 ) + eA 2 (h + 1, k 2 ) (3.15) 
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Using the parametrization (|3.12| ), these equations will give us the energy eigenvalues 

E 2u3+l =Q + Q- 1 + e(6 + & 1 ) (3.16) 

To find £ 2 we must consider the near equations. First, we compute the action of H on the coupled near 
states [Ox (At, fc + 1)) and \Q 2 (k, k + 2)): 

H \fli(k, k + l)) = (Q + Q- 1 ) \Qi(k, k + l))+e \£li(k, k + 2)) + e \Q 2 (k, k + 2)) 

(3.17) 

H |fi 2 (A;, fc + 2)) = (Q + g- 1 ) |fi 2 (fc, k + 2)) + e |ft 2 (A; - 1, A; + 2)) + e |fii(A;, fc + 1)) 

(3.18) 

The last terms in these equations tell us that a pseudoparticle can propagate past the isolated impurity, 
but in so doing causes a shift in its position by two lattice site. Substituting ( |3.17|) and ( |3.18|) into the 
eigenvalue equation, we get 



E 2uJ+1 -Q- Q' 1 ) A x {k, k + 1) = eAi(k, k + 2) + eA 2 (k, k + 2) (3.19) 
Bhu+i -Q- Q' 1 ) A 2 (k, k + 2) = eA 2 (k - 1, k + 2) + eA x {k, k + l) (3.20) 



These equations, which are not automatically satisfied by the ansatz Q3.12|) , are equivalent to the condi- 
tions 

A 1 (k,k) = A 2 (k,k + 2) , A 1 (k,k + l) = A 2 {k + l,k + 2). (3.21) 

obtained by subtracting Eq. ( [3.19|) from Eq.( [3.14|) for k\ = k , k 2 = k + 1 and by subtracting Eq. ( [3.20|) 
from Eq.( |3.15| ) for k\ = k , k 2 = k + 2, respectively. The conditions ( |3.21|) require a modification of the 
amplitude relation (|3.13| ): 

^ = ^ 2 = & =► Z^I = 1 or = 1 (3-22) 

Putting £j = e ldi , i = 1, 2, it means cos(iV — 2)9 2 = cos 2^. Hence 

2irm±4:iTl/N , . 
9 2 = — — , I and m integers. (3.23) 

In other words, ^2^+1(^1,6) are eigenstates of H with energy eigenvalues given by E 2uJ +i = Q + Q^ 1 + 
2ecos6 l 2. Note that when Nm±2l is a multiple of (N — 2) we get states which are degenerate with the 
one-pseudoparticle states ^ 2uj , which lie in the sector r = 2cu. 

In the sectors r = 2cu + I we also will find states, which consist of one pseudoparticle with parameter 
interacting with I impurities, distributing according to (|3.1|), with parameters = 1,2...,/. 

The energy of these states is parametrized as in ( [3.16| ) and satisfies the condition ( |3.22|) with 
£ — 61 " ' ' 6 6+1 ■ It involves only <^ + i and £ imp = £1 £ 2 • • • 6- being therefore highly degenerate, i.e. 

&^-tf = l (3-24) 

This is to be expected due to the irrelevance of the relative distances, up to jumps of two positions via 
exchange with a pseudoparticle. Moreover, these results do not depend on impurity type. 
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3.3 Two-pseudoparticles 

The sector r = Auj contains, in addition to the cases discussed above, states which consist of two 
interacting pseudoparticles. We seek eigenstates in the form 

*4u,(£i,6)= E A^ki) in^h)) (3.25) 

fcl+l<fc 2 

Again, we try to build two-pseudoparticle eigenstates out of translational invariant products of one- 
pseudoparticle excitations at k\ and k 2 (k 2 > k\ +2) : 

^(6, 6) = *2c(6) x + *2«(6) x 'Mk) (3.26) 

Using again (|3~5| ) and comparing ( |3.26|) with (|3.25| ) we get 

a, f3=—u) 

for k 2 > k\ + 3 and 

A(k u k 2 ) = Au&& + A 21 &g*, (3.28) 

for k 2 > fci +2. 

Periodic boundary conditions A(k 2 , N + fci) = fc 2 ) and ^4(iV + k\,N+ k 2 ) = A(ki, k 2 ) imply 

A 12 g = A 21 and £ N = 1 (3.29) 

where £ = £1^2 (6 = e ie % i = 1, 2) and the total momentum is 9i + 9 2 = 2nl/N, with I integer. 

Applying H to the state of Q3.25|) , we obtain a set of equations for the wavefunctions A(ki, k 2 ). When 
the two pseudoparticles are separated, (k 2 > k\ + 3) these are the following far equations: 

(E 4u -2Q-2Q- 1 )A(k 1 ,k 2 ) = eA(h - 1, k 2 ) + eA(h + 1, A*) 

(3.30) 

+ eA(/c 1 ,/c 2 -l) + eA(/c 1 ,fc 2 + l) 

We already know them to be satisfied, if we parametrize A(k\,k 2 ) by plane waves ( [3 .28|) . The corre- 
sponding energy eigenvalue is 

E 4uj = 2Q + 2Q- 1 + e + fr 1 + 6 + £ 2 ~ X ) (3-31) 

The real problem arises of course, when pseudoparticles are neighbors, so that they interact and we 
have no guarantee that the total energy is sum of single pseudoparticle energies. 

Acting of TC on the state ( |3.27|) gives the following set of equations for the near states 

H\n(k,k + 2)) = 2(Q + Q- 1 )\Q(k,k + 2)) + e\n(k-l,k + 2)) 

(3.32) 

+ e\n{k,k + 3)) + U k+1 \Q{k,k + 2)} 

Before we substitute this result into the eigenvalue equation, we observe that some new states are 
appearing. In order to incorporate these new states in the eigenvalue problem, we define 

U k+1 \tt(k, k + 2)) = t \tt(k, k + l)) + e \n(k + l,k + 2)) (3.33) 
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Here we underline that we are using the same notation for these new states. Applying H to them we 
obtain 

n\n(k,k + i)) = (Q + Q-^in^k + i)) + e\n(k-i,k + i)) 

(3.34) 

+ e\Q(k,k + 2)} 
Now, we extend ( |3.25|) , the definition of ^> 4uJ , to 



*4u,(£i,6)= E A(k u k 2 )\n(k u k 2 )} (3.35) 

ki<k 2 



Substituting ( [3.32|) and ( B.34|) into the eigenvalue equation, we obtain the following set of near equations 

(E Au -Q- Q- 1 ) A(k, k + 1) = eA(k - 1, k + 1) + eA(k, k + 2) (3.36) 



Using the same parametrization ( |3.28|) for these new wavefunctions, the equation ( |3.36|) gives us the 
phase shift produced by the interchange of the two interacting pseudoparticles 

A21 e(l+Q + (Q + Q-% 

An ~ e(l+0 + (Q + g-% 1 ■ ' 

We thus arrive to the Bethe Ansatz equations which fix the values of £1 and £ 2 in the energy equation 

^ = = 1 (3.38) 

3.4 Two-pseudoparticles and impurities 

In the sectors r > 6cu, in addition the cases already discussed, we find states with two interacting particles 
and impurities. Let us now consider two pseudoparticles with one impurity of type Theses eigenstates 

lie in the sector r = Auj + 1 and we seek them in the form 

*Wi(&,6,&) = E A 1 (k u k 2 ,h)\n 1 (k u k 2 ,k 3 )) 

/ci+l<fc2<fc 3 -2 

+ E ^2(^1, A; 2 , A; 3 ) |0 2 (/ci, A; 2 , ^ 3 )> 

fci + l<fc 2 <fc3 

+ E Mh,k2,h)\n 3 {k u k2,k 3 )) (3.39) 

fcl+l<fc 2 <fc3-l 

In A{(ki, k 2 , fc 3 ) the index i = 1,2,3 characterizes the impurity position. Comparing ( B.39|) with the 
state build from the translational invariant products of two-pseudoparticles with parameters £ 2 and £3 and 
one-impurity with parameter £x : 

*Wi(&,6,&) = |(w - i)(6)) x *2u,(6) x ^(£3) 

+ |(w-l)(£i))xtt 2w (£g)xtt 2w (£ 2 ) 

+*2 W (6)X|(W-I)(ei)>x* 2w (e3) 
+*2 W (6)X|(W-I)(ei)>x* ai ,(e2) 

+*2,(6)x^fe)x|(a;-l)(a)) 

+* 2w (&) x ^(6) x |(w - 1)(6)) (3.40) 
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we get 



\n 1 (k u k 2 ,k 3 )} = £ W(a,/3,u) \h(uj- l),k 2 (-a,a), k 3 (-/3,P)) 

a, /3=—lj 

\n 2 (h, k 2 , h)} = jr W(a } (3, u) \h(-a, a), k 2 (uj - 1), k 3 (-/3, (3)) 

a, f3=—uj 

\n 3 (k 1 ,k 2 ,h)}= jr wi^Mikxi-^a^hi-^^hiu-i)) (3.4i) 

a, fi=—u) 



where 



^312 


^213 


^231 


^4l32 ' 




= (666) 



W{a, (3, a) = e{a)e{(3) q -< 2 ^+^+^p> (3.42) 
and the wavefunctions Ai(ki, k 2 , k 3 ) which are parametrized by plane waves as 

A 1 (k 1 , k 2 , fc s ) = A 123 6 fcl 6 fc2 6 fc3 + ^i3 2 6 fel 6 fc3 6 fe2 
A 2 {h, k 2 , k 3 ) = A 213 6 fc2 6 fcl 6 fc3 + ^23i6 fc2 6 fc3 6 fcl 

A 3 (h, k 2 , k 3 ) = A^&Zs 2 + ^32i6 fc3 6 fc2 6 fcl - (3-43) 
Periodic boundary conditions read now 

Ai(kx, k 2 , k 3 ) = Ai(N + h,N + k 2 , N + k 3 ), 
Ai(k 2 ,k 3 ,N + h) = A i+1 (k h k 2 ,k 3 ), i= 1,2, 3 mod 3 (3.44) 

which imply that 

_ ^123 _ ^132 £ jV 
^312 ^321 

g = T^ = T^, ^ = (666)^ = 1 (3.45) 

^213 ^123 

Action of H on the state ^4ui+i gives the following set of far equations: 

(e w - 2Q - 2Q- 1 ) A 1 (k 1 , k 2 , k 3 ) = eA 1 {k 1 , k 2 - 1, k 3 ) + eA^kx, k 2 + 1, k 3 ) 

+eA 1 (k 1 , k 2 , k 3 - 1) + eA^h, k 2 , k 3 + 1) 

(3.46) 

and a similar set of eigenvalue equations for A 2 (ki, k 2 , k 3 ) and A 3 {ki, k 2 , k 3 ). The parametrization ( |3.43|) 
solves these far equations provided that 

E w = 2Q + 2Q- 1 + e(6 + + 6 + (3.47) 

Taking into account the near equations we must split them in three different neighborhood: (i) impurity 
neighbors of separated pseudoparticles and, (ii) impurity far from neighbors pseudoparticles and (iii) when 
impurity and pseudoparticles share the same neighborhood. 

In the case (i) we consider the second pseudoparticle far and follow the steps for the case of one- 
pseudoparticle with impurity eigenstates. Thus, the near equations can be read off from ( |3.36| ) 

{E AuJ+l -2Q- 2Q' 1 )A 1 {k, k + 1, k 3 ) = eA 1 {k, k + l,k 3 -l) + eA 1 (k, k + l,k 3 + l) 

+eA 1 (k, k + 2, k 3 ) + eA 2 (k, k + 2, k 3 ) (3.48) 



10 



and a similar set of equations coupling A 2 and A 3 . It follows from the consistency between ( [3.46|) and 
CT ) that 



A 1 (k,k,k 3 )=A 2 (k,k + 2,k 3 ) (3.49) 
and similar identification between A 2 and v4 3 . The plane waves ( |3.43|) solve these identifications provided 



-2 ^123 ^4l32 ^231 ^213 



^ ~~ 4 ~~ 4 ' ~~ 4 ~~ 4 ^ ^ 

^213 ^231 ^321 ^312 

For the case (ii) we can derive the near equations from those of two-pseudoparticles case. Keeping 
the impurity far and following the steps (|3.30| )- (|3.37| ) we get 

Eau+i -Q- Q' 1 ) A 1 (k 1 , k, k + 1) = eA x {kx, k - 1, k + 1) + eA x {k u k,k + 2) (3.51) 



and a similar set of equations for A 2 and A 3 . The case (iii) is obtained from ( |3.51|) for k\ — k — 1. 
The parametrization ( [3.43|) solves this provided that 



^123 e(l + ^ 2 ) + (Q + Q- l )& 

A 132 " e(l + £ 2 £ 3 ) + (Q + q- 1 )^ 1 ' j 

Matching the constraint equations ( |3.52|) , ( |3.50|) and (|3.45|) we arrive to the Bethe equations 

" = -l + « t + E (Q + <?-%' ^ i = 2 ' 3 (3 ' 53) 
£" = «i66)" = 1, ff" 4 = 1. (3.54) 

The origin of the exponent (iV — 4) in the impurity parameter can be understood by saying that after the 
two pseudoparticles propagate past the impurity, the position of impurity is shifted by four lattice sites. 

Next, we can also find eigenstates with two pseudoparticles and more than one impurities. They can 
be described in the following way: Let us consider an eigenstate with I > 1 impurities with parameters 
£i,£2, ■ ■ • >6 and two pseudoparticles with parameters and £/ +2 . The energy eigenvalue is 

E r = 2Q + 2Q- 1 + e + C+i + 6+2 + C+ 2 ) (3-55) 

and the Bethe equations 

i + a+i6+2 + 6(g + g- 1 )6+i 
- &+16+2 + e(g + g- i )a +2 

£f~ 4 = l, a = 1,2,...,/ (3.56) 

Moreover, ^ = 1 with £ = £i£ 2 • • -£z+2- 
3.5 Three-pseudoparticle eigenstates 

In the sector r = 6tu, in addition to the previously discussed eigenstates of one and two pseudoparticles 
with impurities, one can find eigenstates with three interacting pseudoparticles with parameters £1,^2 ar| d 
£ 3 . We start seek them in the form 

&)= E Mki, k 2 , k 3 ) \Q(k l7 k 2 , k 3 )) (3.57) 

fci+2<fc 2 <fe 3 -2 



SZ+1S1S2 * * ' SJ — 
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where k 2 , k 3 )) = ®f =1 \Q(ki)) . The corresponding wavefunctions 

A(k 1 ,k 2 ,k 3l ) = A^e^^+A^^^+A^ 2 ^^ 3 



(3.58) 



+A 231 ei 2 && + A^&zt + A m &&& 

satisfy the periodic boundary conditions 

A(k 2 , k 3 , N + ki) — A(k h k 2 , k 3 ), A(N + k h N + k 2 , N + k 3 ) = A(k 1} k 2 , k 3 ) (3.59) 
which imply that 

An 3 _ A\ 32 ^ N _ A 3 i 2 _ A 2 \ 3 

1 ^312 ^321 ' " ^231 ^132 

c n A 32 \ A 231 N N 

4 3 = ~ A — = ~a — > i = (6424s) = 1 (3.60) 

^213 ^123 

These relations show us that the interchange of two-pseudoparticles is independent of the position of the 
third particle. 

Applying H to (|3.57| ), we obtain a set of equations for A(ki, k 2 , k 3 ). When the three pseudoparticles 
are separated, (kx + 2 < k 2 < k 3 — 2), we get the following far equations: 

(E 6s - 3Q - ZQ-^Afa, k 2 , k 3 ) = eA{k 1 -l,k2,k 3 )+eA{k 1 + l,k 2 ,k3) 

+eA(k h k 2 - 1, k 3 ) + eA(kx, k 2 + 1, k 3 ) 
+eA(k u k 2 , k 3 -l) + eA(kx, k 2 , k 3 + 1) 

(3.61) 

It is simple verify that the wavefunctions (|3.58| ) satisfy these far equations provided 

= £ {Q + Q- 1 + £ n + C 1 } (3.62) 

n=l 

Applying H on the near states we get the following set equations: 

H\n(k 1 ,k 1 + 2,k 3 )) = (2Q + 2Q- 1 )\n(k 1 ,k 1 + 2,k 3 )) + e\n(k 1 -l,k 1 + 2 J k 3 )} 

+e \n{kx, h + 3, k 3 )} + e |0(/ci, k x + 2,k 3 - 1)) 
+e \n(h, h + 2,k 3 + 1)) + U kl+1 \Q{h, h + 2, fc 3 )> 

(3.63) 

for k 3 > ki + 4, which correspond to the meeting of two pseudoparticles at the left of the third pseu- 
doparticle, which is far from of the meeting position. 

HMk 1 ,k 2 ,k 2 + 2)) = (2Q + 2Q- 1 ) \Q{k 1 ,k 2 ,k 2 + 2)) + e^kx - l,k 2 ,k 2 + 2)) 

+e \Q(ki + 1, k 2 , k 2 + 2)) + e \Q(ki, k 2 -l,k 2 + 2)) 
+e \Q(kx, k 2 , k 2 + 3)) + U k2+1 k 2 , k 2 + 2)) 

(3.64) 
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for k 2 > ki + 2, which correspond to the meeting of two pseudoparticles at the right of the far pseudopar- 
ticle. Moreover, there is one set of equations which correspond to the meeting of three pseudoparticles 

H\n(k,k + 2 } k + 4:)} = (Q + Q- 1 )\Vl(k,k + 2 } k + 4:)} + e\n(k-l,k + 2,k + 4)} 

+e \Q(k, k + 2,k + 5)} + U k+1 \Q(k, k + 2, k + 4)) 
+U k+3 \Q(k,k + 2,k + A)) (3.65) 

In deriving these equations new states made their debut. In order to incorporate these new states in the 
eigenvalue problem we define: 

t^ fcl+1 10(^,^ + 2,^)) = e\n(k 1 ,k 1 + l,k 3 )) + e\n(k 1 + l,k 1 + 2,k 3 )) 

U k2+1 \Q(k h k 2 ,k 2 + 2)) = e\n(k 1 ,k 2 + l,k 2 + 2)} + e\n(k 1 ,k 2 ,k 2 + l)} 

U k+1 \Q(k,k + 2,k + i)) = e \Q(k, k + 1, k + 4)) + e \Q(k + 1, k + 2, jfe + 4)) 

U k+3 \n(k,k + 2,k + A)) = e\n(k,k + 3,k + A)) + e\n(k,k + 2,k + 3)) (3.66) 

Applying H to these new states the result can be incorporated to the eigenvalue problem provided the 
definition of ( |3.57|) is extended to 

<M6,6,6) = E Mkuh^Mh,^)) (3.67) 

ki<k 2 <k 3 

After this we are left with three meeting equations 

(E 6uj -2Q- 2Q- 1 )A(k u h + 1, k 3 ) = eA(h -l,h + 1, k 3 ) + eA(k u k x + 2, k 3 ) 

+ eA(k u h + 1, k 3 - 1) + eA(k u h + l,k 3 + 1) 



(3.68) 



for k 3 > ki + 2, 

(E 6uj -2Q- 2Q~ 1 )A(k 1 , k 2 , k 2 + l)= eA(k u k 2 -l,k 2 + l) + eA(k h k 2 , k 2 + 2) 

+ eA(ki - 1, k 2 , k 2 + 1) + eA{k x + 1, fc 2 , fc 2 + 1) 

for ki + 2 < k 2 and 

(E 6u} -Q- Q~ 1 )A(k, k + 1, k + 2) = eA(A; - 1, jfe + 1, jfc + 2) + eA(k, k + 1, jfe + 3) 
It is easy to verify that the parametrization ( |3.58|) and (|3.62| ) solve these equations provided 



(3.69) 



A\ 23 


^231 


i + 66 + 


e(Q + Q" 


- x )6 


^213 


^321 


i + 66 + 


e(Q + Q" 


- x )6 


^132 


^213 


1 + 66 + 


e(Q + Q- 


- x )6 


^231 


^312 


1 + 66 + 


e{Q + Q~ 


- x )6 


^312 


^4l23 


1 + 66 + 


e(Q + Q~ 


"% 


^321 


^4l32 


i + 66 + 


e(Q + Q- 


- x )6 



(3.70) 

Matching these constraints and the periodic boundary conditions fl3.60|) we get the Bethe Ansatz equations 



i + 66 + e(Q + Q" 1 )6 



(666r = i (3.71) 
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3.6 General eigenstates 



The generalization to any r is now immediate. Since the Yang-Baxter equations are satisfied, there is 
only two-pseudoparticle scattering (using the S'-matrix language). Therefore, neighbor equations, where 
more then two pseudoparticles become neighbors, are nor expected to give any new restrictions. For 
instance, in the sector r = 6lu, we saw that the interchange of two-pseudoparticles is independent of the 
position of the third particle. Thus, in a sector with p pseudoparticles we expect that the p-pseudoparticle 
phase shift will be a sum of p(p — l)/2 two-pseudoparticle phase shift. The energy is given by the 
sum of single pseudoparticle energies. The corresponding Bethe Ansatz equations depend on the phase 
shift of two pseudoparticles and on the number of impurity. For a generic sector one can verify that no 
different neighborhood those discussed above can appear. So, no additional meeting conditions will be 
encountered. Thus, we can extend the previous results to the p -pseudoparticle states in the following way: 
In a generic sector r with I impurities parametrized by £162 •• "6 and p pseudoparticles with parameters 
6+16+2 ••• £i+p, the energy is 

E r = E {Q + Q^ + e^n + C 1 )} ( 3 - 72 ) 

n=l+l 

with £ n determined by the Bethe ansatz equations 

fNf 2 f 2 ,2 yf f i + ^a + e(g + g-% l 

{aClCa 1 " J-i+il l + U b + e(Q + Q-^ b f 

ic~ 2v = 1, c=l,2,...,Z 

i N = 1, £ = £i£ 2 -- -66+16+2 •••&+*, (3.73) 

The energy eigenvalues and the Bethe equations depend on the deformation parameter q, through the 
relation ( PTT|) : 

r [2s + 1] for Ai 

n + n-i- [2n-l][n+l/2]/[n-l/2] for B n (n > 2) 

^ + ^ I [n][2n + 2]/[n + l] for C n (n > 1) ^' J 

[ [n)[2n - 2)/[n - 1] for D n (n > 3) 

We obtained thus the spectra with periodic boundary conditions of quantum spin-chain models, arising 
as representations of the Temperley-Lieb algebra. As expected, all these models have equivalent spectra 
up to degeneracies of their eigenvalues. From a suitable sorting of the parameters £ it one can insure that 
the spectra of lower-r sectors are contained entirely in the higher-r sectors. 



4 Bethe Ansatz: Non-local boundary conditions 



It is the purpose of this section to present and solve, via coordinate Bethe ansatz, the quantum group 



invariant closed TL Hamiltonians which can be written as 19 



N-l 

H= E Uk+Uo 
fe=i 



(4.1) 
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where U k is a Temperley-Lieb operator and bt = U is non-local term defined through of a operator G 
which plays the role of the translation operator 

U = GU N ^G~ l , G = (Q-U 1 )(Q-U 2 )---(Q-U N - 1 ) (4.2) 

satisfying [H,G] = and additionally invariance with respect to the quantum algebra. The operator 
G shifts the U k by one unit GUkG^ 1 = U k +i and maps Uq into U\, which manifest the translational 
invariance of H. In this sense the Hamiltonian (|4~1] ) is periodic. 

For the g-deformed A-D Temperley-Lieb algebra, the matrix elements of U k is again given by (|2~8|), 
i.e. : 

(u,X\U\fi, k) = e(//)e(v)r <e * +e "' p> W<WA,o 

k, A E J (4.3) 

From ( |4.3|) we choose an particular Bethe state 

|0(*)> = £ e (w)e(a) g~< e ^> |*(-a,a)> (4.4) 

a=— lu 

which is an eigenstate of Uk and it is shifted by one unit under the action of Uf.±i 



U k \{l(k)) = (Q + Q- 1 )\tt(k)) 

u k±1 \n(k)) = e \n(k ± i)) , u k \n(k±)) = e \n(k)) 

U k =0 for k {I - 1, I, I + 1} (4.5) 

where e = 1 for B n , D n and Ai ( s integer) and e = — 1 for C n and Ai ( s semi-integer). 

The action of the operator G on the states | can be easily computed using (|4~5| ): It is simple on 
the bulk and at the left boundary 

G \Q{k)) = -eQ N - 2 \Q(k + 1)) , l<£;<iV-2 (4.6) 
but manifests its nonlocality at the right boundary 

JV-l 

G \n(N - 1)) = eQ N ~ 2 J2 (~tQY k P(N - k)} (4.7) 

k=\ 

Similarly, the action of the operator G^ 1 = (Q^ 1 — U^-i) • • • (Q^ 1 — U\) is simple on the bulk and at 
the right boundary 

G- 1 \Q(k)) = -eQ- N+2 \Q(k - 1)) , 2<k<N-l (4.8) 
and non-local at the left boundary 

N-l 

G- 1 |n(l)> = eQ- N+2 !«(*)> ■ (4-9) 

k=l 

Now we proceed the diagonalization of H as was made for the periodic case. As and ( p..2|) have 
the same bulk, i.e., differences came from the boundary terms, we will keep all results relating to the bulk 
of the periodic case presented in the previous section. 
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4.1 One-pseudoparticle eigenstates 

Let us consider one free pseudoparticle which lies in the sector r = 2uj 



N-l 

*2u> = E A (k) Mk)) . (4.10) 
k=l 

The action of the operator U = Y,k=i £4 on the states \Q(k)) can be computed using ( |4.5|) : 

U\Q(1)) = (Q + Q- 1 ) |fl(l)) + e|n(2)) 
U \Q{k)) ={Q + Q- 1 ) \Q(k)) + e \Q(k - 1)) + e \Q(k + 1)) 

for 2 < fc < iV - 2 

W |0(AT - 1)) = (Q + Q- 1 ) \n(N - 1)) + e |fi(JV - 2)) . (4.11) 
and using (|4.6|)-(fO|) one can see that the action of Wo = GUn-iG^ 1 vanishes on the bulk 

U \tt(k)} = , 2<k<N-2 (4.12) 

and is nonlocal at the boundaries 

N-l N-l 

U |0(1)> = -e £ (-eQ) fe |fi(fc)> , U \Sl(N - 1)) = -e £ (-eQ)-" + * |fi(fc)> ■ (4.13) 

k=l k=l 

which are connected by 

Uo \Q(N - 1)) = (-eQ)- N U |fi(l)> • (4.14) 

From these equations we can understood the role of U : Although the Hamiltonian (|4.1|) is a global 
operator, it manifests the property of essential locality. From the physical point of view, this type of 
models exhibit behavior similar to closed chains with twisted boundary conditions. 

Before we substitute these results into the eigenvalue equation, we will define two new states 

e |fi(0)) = U Ml)} , e \n(N)) = U \Q(N - 1)) (4.15) 



to include the cases k = and k = N into the definition of equation ( |4.10|) . Finally, the action of 
H = U + Uq on the states \Vt(k)) is 



H |fi(0)) = (Q + Q- 1 ) MO)) + (-eQ) iV e \Sl(N - 1)} + e 

W TO) = (Q + g- 1 ) |fi(A;)) + e - 1)> + e + 1)) 

for 1 < jfe < JV - 2 

W |fi(JV - 1)) = (Q + g^ 1 ) |fi(JV - 1)) + e |fi(JV - 2)) + (-eQ)~ N e |fi(0)) 

W |n(JV)> = (Q + g _1 ) |fi(A0) + e |0(AT - 1)) + (-eQ)- JV e (4.16) 
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Substituting these results into the eigenvalue equation H^ 2uJ = E 2u} ^ 2ui we g et a complete set of 
eigenvalue equations for the wavefunctions 

E 2s A(k) = (Q + Q- l )A(k) + eA(k - 1) + eA(k + 1) 

for 1 < k < N - 1 (4.17) 

provided the following boundary conditions 

{-eQ) N A{k) = A(N + k) (4.18) 

are satisfied. 

The plane wave parametrization A(k) = A£ k solves these eigenvalue equations and the boundary 
conditions provided that: 

E 2LU = Q + Q- l + e{i + C 1 ) 

e = (-eQf (4.19) 

where £ = e ld and 9 being the momentum. 
4.2 Two-pseudoparticle eigenstates 

Let us now consider the sector r = Auj, where we can find an eigenstate with two interacting pseudopar- 
ticles. We seek the corresponding eigenfunction as products of single pseudoparticles eigenfunctions, 
i.e. 

*4u= E Mk\,kz) \Q(k 1 ,k 2 )) (4.20) 

k!+Kk 2 

where 

|n(*i, k 2 )) = jr e(a)e((3) q -<^ + e a+ep , P > a)> p)) (4 21) 

a, [3=—lu 

To solve the eigenvalue equation = E^^^, we recall ( |4.5| ) to get the action of U and Uq 

on the states \VL{k\,k 2 )) . Here we have to consider four cases: (i) when the two pseudoparticles are 
separated in the bulk, the action of U is 

U \Q(h, k 2 )) = 2(Q + Q- 1 ) k 2 )) + e \n(k x - 1, k 2 )) + e \Q(h + 1, k 2 )) 

+e \Q(ki, k 2 - 1)> + e \Q{k u k 2 + 1)) (4.22) 

i.e., for ki > 2 and ki + 3 < k 2 < N — 2; (ii) when the two pseudoparticles are separated but one of 
them or both are at the boundaries 

U fc 2 )) = 2{Q + g- 1 ) |0(1, k 2 )) + e |fi(2, fc 2 )) + e A; 2 - 1)) 

+e |0(1, A; 2 + 1)) (4.23) 

U |0(A;i, iV - 1)) = 2(Q + Q" 1 ) |0(A;i, N - 1)) + e ^(^ - 1, JV - 1)) 

+e + 1, iV — 1)) + e ^(jfei, iV - 2)) (4.24) 
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U N — i)) — 2(Q + Q- 1 ) iV — 1)) + e |H(2, TV - 1)) + e iV - 2)) (4.25) 
where 2 < fci < iV — 4 and 4 < k 2 < N — 2; (iii) when the two pseudoparticles are neighbors in the bulk 

U \n(k, k + 2)) — 2(Q + Q' 1 ) \Q(k, k + 2)) + e \Q(k -l,k + 2)) + e \Q(k, k + 3)) 

+U k+1 \Q(k,k + 2)) (4.26) 

for 2 < k < N — 4 and (iv) when the two pseudoparticles are neighbors and at the boundaries 

U 3)) = 2(Q + Q" 1 ) 3)) + e |0(1, 4)) + C/ 2 |Q(1, 3)) (4.27) 

W |H(iV - 3, N - 1)) = 2(Q + Q" 1 ) |ft(7V - 3, TV - 1)) + e |f2(iV - 4, TV - 1)) 

+C/ JV _ 2 |fi(iV-3,JV-l)) (4.28) 

Moreover, the action of U does not depend on the pseudoparticles are neither separated nor neighbors. 
It is vanishes in the bulk 

U \Q(ki, k 2 )) = for h ^ 1 &ndk 2 ^ N - 1, (4.29) 
and different of zero at the boundaries: 

U \Q(l,k 2 )) = - e 1] 2 (-eQ) fc |^(^^))-(-eg) fc2 ~ 1 ^ 2 |fi(A;2-l,A; 2 +l)) 

k=i 

N-l 

-e £ (-eQ) k ~ 2 m 2 ,k)} (4.30) 

k=k 2 +2 

U N-l)) = {-eQT N+2 U fc 2 )) (4.31) 

where 2 < fc x < N - 3 and 3 < k 2 < N - 2. 

Following the same procedure of one-pseudoparticle case we again define new states in order to have 
consistency between bulk and boundaries terms 

UoMl,k 2 )) = e|ft(0,fc 2 )), U \n(k 1 ,N-l)) = e\n(k 1 ,N)) 
U \n(l,N-l)} = e\n(0,N-l)) + e\n(l,N)) 
U k+1 \Q(k,k + 2)) = e\n(k,k + l)) + e\£l(k + l,k + 2)) (4.32) 

Acting with H on these new states, we get 

H\n(0,k 2 )) = 2(Q + Q' 1 )\Q(0,k 2 )) + e\Q(0,k 2 -l)) + e\n(0,k 2 + l)) 

+e k 2 )) + (-eQf- 2 e \Q(k 2 , N - 1)> (4.33) 

H \Q(ki, N)) = 2(Q + Q- 1 )\n(k 1 ,N)) + e ^(A* - l,iV)) + e^h + 1, N)) 

+e \Q(h, N-l)) + (- € Q)- N+2 e |0(1, fcx)> (4.34) 

W |ft(Jfe, fc + 1) = (Q + Q' 1 ) \n(k, k + 1) + e |n(A; - 1, A; + 1) + e |ft(Jfe, k + 2) (4.35) 
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Substituting these results into the eigenvalue equation, we get the following equations for wavefunctions 
corresponding to the separated pseudoparticles. 

{E 4w -2Q -2Q- 1 )A(k 1 ,k 2 ) = eA(k 1 -l,k 2 )+eA(k 1 + l,k 2 ) 

+eA(k l ,k 2 -l) + eA(k 1 ,k 2 + l) (4.36) 

i.e., for ki > 1 and k\ + 3 < k 2 < N — 1. The boundary conditions read now 

A(k 2 , N + ki) = (-eQ) N - 2 A(k l7 k 2 ). (4.37) 

The parametrization for the wavefunctions 

A(k u k 2 ) = A l2 ^i k 2 2 + Mx&& (4.38) 
solves the equation (|4.36| ) provided that 

E is = 2{Q + Q~ l ) + e(6 + + 6 + (4-39) 



and the boundary conditions ( |4.37|) provided that 



& = i-eQf" 2 ^ , & = (-eQf- 2 ^ =► e = (-eQ) 2 ^ (4.40) 

A 12 A 21 

where £ = £i£ 2 = e i{8l+d2) , ^ + 6» 2 being the total momenta. 

Now we include the new states ( |4.32j ) into the definition of in order to extend (|4.20| ) to 

*4c= E Aiku^Mh,^). (4.41) 

fcl<fc2 

Here we have used the same notation for separated and neighboring states. 
Substituting ( |4.26|) and (|4.35| ) into the eigenvalue equation, we get 



(E4,, -Q- Q-^Aik, k + l) = eA{k -l,k + l) + eA{k, k + 2) (4.42) 
which gives us the phase shift produced by the interchange of the two pseudoparticles 

A 2 i i+e + e(Q + Q _1 )6 



A 12 i+f + e (Q + 3-1)6- 
We thus arrive to the Bethe ansatz equations which fix the values of £x and £ 2 : 



(4.43) 



{ l + ^ + e(Q + g- 1 )6j 
= (-eQ) 2(Af ~ 2) (4.44) 
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4.3 General eigenstates 

Thus in the sector r = 2ujp, we expect that the p-pseudoparticle phase shift will be a sum of two- 
pseudoparticle phase shifts and the energy is given by 

E p{2s) = Y J {Q + Q- 1 + e{i n + C 1 )} (4-45) 

n=l 

where 

= i-^Q) K -* a n (- ! t ff t In t 1 ■ • = 1 --* 

(ei6---e P ) JV = (-egr (JV - 2p+2) (4.46) 

The corresponding eigenstates are 

*r(£l, 6, -4) = E ^ fc 2,..., fcp) fc 2, fcp)> (4.47) 

l<fei<...<fcp<AT-l 

where k%, = ®f =1 [O(fcf)) and the wavefunctions satisfy the following boundary conditions 

fc 2 ,..., fcp, iV + fci) = (-eQ) JV - 2p+2 A(A; 1 , fc 2 ,.., fc p ) (4.48) 

It is not all, in a sector r we may have p pseudoparticle and N^-i, iV w _2, -W-uj+i impurities of the 
type (a> — 1), (a; — 2), (— + 1), respectively, such that 

AUi + 2A^ 2 + • • • + (2cu - l)N^+i = r-2up (4.49) 

We called impurity a state \a, k) flanked by at least two states |/3, k ± 1) such that a + /3 7^ 0. Since 7i 
is a sum of projectors on spin zero, these states are also annihilated by U . Therefore the impurities play 
here the same role as in the periodic case. It means that for a sector r with I impurities with parameters 
£i,...,£z and p pseudoparticles with parameters £z+i, £z+ p the energy is given by ( ft.46| ), and the Bethe 
equations do not depend on impurity type and are given by 

^Zti+x I 1 + &»6 + e(g + Q l )i b J 

with a = I + 1, f* + 2, f* + p , p > 1, and 

£ 2p & + i • • • 6 + X~ 2p = (-eQ)*"-^ (4.51) 

where £ = fi6j--*66+i 

We have shown that these closed Temperley-Lieb quantum invariant spin chains can be solved by the 
coordinate Bethe ansatz. A consequence of the nonlocal terms U is the arising of boundary conditions 
depending on the quantum group parameter q via the relation Q + Q^ 1 =Try A (q~ 2p ). It is also p- 
pseudoparticle dependent (which is equal to spin sector r for A\, when s = 1/2). 

For the algebra A\ } and s = 1/2 , Q=q and £4 are 4x4 matrices giving a nearest-neighbour interaction 
Uk = q-alo k+1 -o k o+ +1 -(q + q- 1 )/A(a z k a z k+1 + l) + (q + q- 1 )/A(a z k -a z k+l -2). This Hamiltonian was 
investigated in Ref . pO] and we summarize some basic results regarding this case. Assuming q=exp(i</?) 
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some interesting properties were found. For instance, the spin L of the ground state becomes ip - 
dependent. For any N (even), L depends on the value of ip according to : 



L = for — < ip < 7i 

L = 1 for 2ffTl) < V < 27 

N r 7T 

L = — for < w < — 

2 iV 

The ground state is non-degenerate (up to the trivial SU 9 degeneracy). At the edges of the intervals, 
cp = tt/21, additional degeneracies occur. These transitions to higher spins resemble the incommensurate 
transition obtained in various other models. 

From the statistical mechanics point of view the Hamiltonian presents critical behavior and it is 
conformal invariant. The central charge ( or conformal anomaly ) is : 

Tiip 

In particular, if we choose the rational form : 

Tim 

ip = - — , m = 3,4..., 

[m + 1) 

then 

C = 1 ; — , 



777.1771 



which give us the conformal anomalies of the minimal unitary models. 

The connection between the Hamiltonian of the closed SU q (2) invariant chain and the unitary minimal 
series was explored in [f35[]. For a generic irrational cp one can decompose the space of states into the 
direct sum of irreducible representations of the quantum group which are in one-to-one correspondence 
with the usual SU(2) representations. 

We have constructed and diagonalized numerically the Hamiltonian for small values of N and s=l/2, 
1. We checked that, for a given N , s=l/2 and 1 have the same spectra up to degeneracies. 



5 Bethe Ansatz: Free boundary conditions 



It is for free boundary conditions that the Hamiltonian (|2.7| ) naturally commutes with the quantum group 
U q (X n ). Since the our linear combination ( ^2| ) left all models with the same status, which concern to 
the coordinate Bethe ansatz, we expect that all procedure developed for the coordinate Bethe ansatz with 
free boundary conditions in || for the case A\ (s = 1/2). can be used here. To show this we recall 
the previous section, taking into account Uq = 0, where almost all equations can be seized for the free 
boundary conditions eigenvalue problem. 
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5.1 One-pseudoparticle 

In this sector, the eigenstate is given by ( |4.10|) : 

jV-1 

<mo = E m \m) (5.i) 

k=l 

where \£l(k)) is again given by 

The action of H on the states \Q(k)) is given by ( |4.11| ), which gives us the following eigenvalue 
equations 

(E 2u -Q-Q- 1 )A(k)=eA(k-l) + eA(k + l), 2 < k < N - 2 (5.2) 
At the boundaries, we get more two slightly different equations 

-Q- Q-^Ail) = eA(2) 
{E 2ui - Q - Q- l )A{N - 1) = eA(N - 2) (5.3) 

where e = 1 for B n , D n and A\ ( s integer) and e = — 1 for C n and A\ ( s semi-integer). We now try as 
a solution 

A{k) = K(6)i k - A(-8)C k (5.4) 

where £ = e l6 ' , 9 being the momenta. Substituting this in equation ( |5.2|) we obtain the energy eigenvalue 
associated with a free pseudoparticle with free boundary conditions 

^ = Q + Q~ 1 + e(e + r 1 ) (5-5) 

We want equations ( |5.2| ) to be valid for k = 1 and k = N — 1 also, where A(0) and A(N) are defined 
by (|574|). Matching ( |5.2|) and ( |573| ) we get the end conditions 

A(0) = and A(N) = (5.6) 

implying that A{6) = A(— 9) and ^ 2N = 1, respectively. A{6) it now determined ( up to a factor that is 
invariant under 9 < — > —9), to be equal to !;~ N . 

5.2 One pseudoparticle and impurities 

Differently from the previous cases, due to the lack of periodicity, the impurity positions are fixed. So, they 
have a different role in the eigenvalue problem with free boundary conditions. For instance, let us consider 
the case of one impurity of the type u — 1, with parameter £x and one pseudoparticle with parameter £ 2 - 
This eigenstate lies in the sector r = 2u + 1 and we can write 

* 2w+ i(ei,6)= E {A 1 (k 1 ,k 2 )\n 1 (k u k 2 )) + A 2 (k u k 2 )\n 2 (k u k 2 )}} (5.7) 

k±<k2 

where |fii(/ci, k 2 )) , % — 1, 2 are given by (|3.11| ). 

For this case we obtain the following eigenvalue equations 

(E 2uj +i - Q-Q~ 1 )A 1 (k 1 ,k 2 ) = eA 1 (k 1 -l,k 2 ) + eA 1 (k 1 + l,k 2 ) 
{E 2bJ+x -Q-Q- l )A 2 {k x ,k 2 ) = eA 2 {k 1 ,k 2 -l) J reA 2 {k l ,k 2 + ^-) (5.8) 
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We also have two meeting conditions that arise because pseudoparticle and impurity may be neighbors 
(see ([J]?!])) 

Ax(k, k) = A 2 (k, k + 2), A 2 (k + l,k + 2) = A t (k, k + 1) (5.9) 
in addition to the two conditions to be satisfied at the free ends 

A 1 {k 1 ,N) = , A 2 {0,k 2 )=0 (5.10) 

Now we try the following ansatz for the wavefunctions 

A 1 (k 1 ,k 2 ) = A 1 (9 1 ,9 2 )^-A 1 (9 1 ,-9 2 )^ 2 ^ 

A 2 (h,k 2 ) = A 2 {9x,9 2 )^^-A 2 {9 u -9 2 )i 2 k ^^ (5.11) 
From ( |5.8|) we get the energy eigenvalue 

£Wi = Q + Q' 1 + e (6 + (5.12) 

and from ( |5.9|) and ( [5.1 0|) the following relations between the coefficients Ai 

A 1 (9 l ,9 2 )^ = A 1 (9 1 ,-9 2 )£ 2 N , A 2 (9 1 ,9 2 )=A 2 (9 1 ,-9 2 ) 

A 1 {9 1 ,9 2 ) = A 2 (9 ll 9 2 )^l , A 1 (9 1 ,-9 2 )=A 2 (9 1 ,-9 2 )e i (5.13) 

from this we get 

e 2 N = i (5.i4) 

as the Bethe equation of (|5.12j ). The coefficients Aj are determined up to a factor that is invariant under 



9 2 < — > —9 2 as: 

Ai(0i, 9 2 ) = e^2 N and A 2 (9 1 , 9 2 ) = & N . (5.15) 

In general, for the eigenstate with I impurities with parameters £i, and one pseudoparticle with 
parameter which lies in a sector r, we can write 

*r(£i, = ±\ E fc+i) 1 (5-16) 

i=l [l<fei<...<fe i+1 <Af-l J 

The corresponding eigenvalue is given by ( |5.5|) , with £ = and the ansatz for the coefficients of the 
wavefunctions becomes 

^u--'A + i)=yu^^H ( 5 - 17 ) 

Here we notice that the index j in the wavefunctions Aj(ki, ki+i) means that the pseudoparticle is at 
the position h +2 -j. 

5.3 Two-pseudoparticles 

For the sector r = 4u, beside eigenstates with impurities, we have an eigenstate with two pseudoparticles. 
We obtain the following eigenvalue equations 

(E 4uj -2Q- 2Q~ 1 )A(k 1 , k 2 ) = eA(h - 1, k 2 ) + eA(h + 1, k 2 ) 

+eA(k l ,k 2 -l) + eA(k 1 ,k 2 + l) (5.18) 
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We have again two conditions to be satisfied at the ends of the chain 

A(0,k 2 )=0 and A(h, N) = (5.19) 
In addition to this we have a meeting condition 

eA(k, k) + eA(k + l,k + l) + (Q + Q~ 1 )A{k 1 k + 1) = (5.20) 
Now we try the ansatz 



A(k u k 2 ) = A(e 1 ,d 2 )^2 2 -Ho 2 ,e 1 )^ 

A(-e u e 2 )^& + A(-e 2 , e,)^^ 
a(9 1 , -9 2 )ei^2 k2 + A(e 2 , -e^ 2 ^ 1 
+ A{-e u -«r fel 6~ fc2 - A(-e 2 , -#i)6~ fc2 6~ fel 



(5.21) 



Here we observe the permutations and negations of 9i and 9 2 . Substituting this ansatz in ( |5.18|) we obtain 
the energy eigenvalue for the sector with two pseudoparticles 



E Auj = 2Q + 2Q- 1 + e (Ci + 6" 1 + 6 + (5.22) 

The ansatz (|5.21j ) satisfy equations ( |5.19|) provided that 

A(9 1 ,9 2 )=A(-9 1 ,9 2 ) , A(9 2 ,9 1 )=A(-9 2 ,9 1 ) 

A(6 1 , -9 2 ) = A{-9 U -9 2 ) , A(9 2 , -9,) = A(-9 2 , -9,) (5.23) 



and 



A(9 l} -9 2 ) A(-9 l} -9 2 ) N A(9 2 ,-9 1 ) A(-9 2} -9,) 
^ 2 " A(9 U 9 2 ) A(-9 U 9 2 ) ' ^ " A(9 2l 9 l ) A{-0 2 ,9 X ) 1 ' ) 

Moreover, the meeting conditions are satisfied provided that 

A(-0i, -O2) _ A(9 2 , gi) = 1 + 66 + 6 (Q + Q- 1 ) 6 
A(-0 2 , -0 X ) A(9 U 9 2 ) 1 + 66 + e (Q + Q- 1 ) 6 
A(9 1 ,-9 2 ) _ A(9 2 ,-9 1 ) l + 6" 1 6 + e(g + g- 1 )6 r , 2 ,x 

a(-9 2 ,9 1 ) A(-e 1 ,e 2 ) i + 6 _1 6 + e(Q + g- 1 )6" 1 

Matching these conditions we get 

2N _ B(-6 1 ,6 2 ) 2N _ B(-8 2 ,9i) , , 

4l ~ B(9 U 9 2 ) ' 42 " B(9 2 ,9 X ) 

and 

A(9 U 9 2 ) = ^ N ^ N B(-9 U 9 2 )&. (5.27) 
Here we have used the usual free boundary notations 

B(9 a ,9 b ) = s(9 a ,9 b ) s(9 b ,-9 a ) (5.28) 

where 

s(9 a , 9 b ) = l + 66 + e (Q + Q' 1 ) 6- (5.29) 
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Now let us consider the eigenstates with two pseudoparticle and impurities. The energy eigenvalue is 
the same of the two pseudoparticles pure state. The parameters associated with impurities are embraced 
in the definition of the coefficients of the wavef unctions. For instance, when we have an eigenstate of two 
pseudoparticles with parameters £ 2 and £3 and one impurity of parameter £ 1( the energy is given by ( (b.22| ) 
and the Bethe equations by ( |5.26| ), with £ x — > £ 3 and Q\ — > #3. But now the wavefunctions are different 

Ai (0i A, 3 ) = {t$t 2 N t 3 N B{-9 2 ,9 3 )& 
A 2 (6 l ,e 2 ,e 3 ) = {^^^Bi-e^e,)^ 1 

A 3 (0i, 2 , 9 3 ) = & N t 3 - N B{-6 2 , e 3 )^ (5.30) 
where B(-9 2 ,9 3 ) is given by (|]28|) 



5.4 General eigenstates 



The generalization follows as in the previous cases. In a sector r with 
p pseudoparticles, we get 



n=l 



and the Bethe equations 



2N 



l+p 

n 



B(—9 a , 9b 



The corresponding eigenfunction can be written as 



a = 1,2,..., 



V 



fci<---<fe i+p 



h +p ) |fi(fci, k 2 , k p )) 



with 



A(h, k 2 , ...,k p ) = E e P A(0 1; e 2 , e p ) ■ ■ ■ # 



(5.31) 



(5.32) 



(5.33) 



(5.34) 



where the sum extends over all permutations and negations of 9 1 , ...,9 P and ep changes sign at each such 
interchange. The coefficients in the wavefunctions are given by 



A(0i, 2 , 9 P 



1R v II s( 9. r o^ x 

j=l l+l<j<i<l+p 



(5.35) 



where B(-9 J ,9 l ) are defined in ( ^28|) . 

For a sector r with I impurities with parameters and p pseudoparticles with parameters 

...,^+p the energy is given by ( |5.31|) and the Bethe equations by ( |5.32|) . Only the coefficients of the 
wave functions are modified 

Aj(0i, 2 , 9 l+p ) = Ajfafr ■ ■ ■ 6)A(0 i+ i, 2 , 0* +p ). (5.36) 

The functions Aj(£i£ 2 ■ ■ -^) = ^i 1 ^ 2 ' ' where the index j characterizes the possible configurations 
of I impurities relative to the p pseudoparticles. Here are numbers which depend on the position of 
corresponding impurity relative to the pseudoparticles. 

Here we observe again the valid of these results for all Temperley-Lieb spin chain Hamiltonians defined 
as projector of spin zero on the representations of the quantum groups U q (X n ), characterized by the values 
of Q + Q- 1 = Ti VA (q- 2p ). 
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6 Conclusion 



We have applied in a systematic way the coordinate Bethe ansatz to find the spectra of a series of "spin 
Hamiltonians arising as representations of the Temperley-Lieb algebra. We consider several boundary 
conditions in order to include all previously known cases. 

Due to Uk be a projector of spin zero, there is a linear combination of eigenstates of Sf = Y,k Sf 
where Sf = diag(a;, — u), to = max(J), which beside simplify the calculus permits us a unified 
treatment for all models. We find that for a given set of boundary conditions, all models have equivalent 
spectra, i.e. they differ at most in their degeneracies. Moreover, for the closed cases, the spectra of the 
lower-dimensional representations are entirely contained in the higher-dimensional ones (see Eq.( |3.23| )). 

Here we notice that this spectrum equivalence is, of course, a consequence of the TL algebra. Never- 
theless there is in the literature a large class of Hamiltonians which are not derived from representations 
of the TL algebra which share the same property. The authors of reference developed a technique for 
construction of spin chain Hamiltonians which affine quantum group symmetry whose spectra coincides 
with the spectra of spin chain Hamiltonians which have non-affine quantum group symmetry. 

The energy eigenvalues are given by 



E = Y,(Q + Q- 1 + 2ecos6 r 



n=l 



Q + Q- 1 = Tr yA (g- 2 (6.1) 

where p is half the sum of the positive roots of U q (X n ), X n = Ai, B n , C n and D n . n are solutions of 
the Bethe ansatz equations (|373]), ( |43Q|) and (ggg ). 



The Hamiltonians for the cases X n = B n , C n and D n appear to be new although due to the Temperley- 
Lieb equivalence 0], they are expected to possess the same thermodynamic properties as the A\(s = 1/2) 
case, i.e., the spin-1/2 XXZ chain with appropriate coupling. 

There are several issues left for future work. In particular, one would like clarify from an algebraic 
point of view, the equality of the spectra, for instance, of the biquadratic model (Ai,s = 1) and XXZ 
model (A 1 ,s = 1/2 ) for free boundary conditions and the inclusion of the XXZ spectrum in the one 
of the biquadratic Hamiltonian for periodic boundary conditions. Furthermore, the completeness and 
complete characterization as highest weight states of the Bethe ansatz eigenstates here presented are not 
considered. 

Using our solutions, one can derive partition functions in the finite-size scaling limit and find the 
operator content of the systems constructed from these quantum chains. 

Finally, we remark here that although the Hamiltonian ( ft.lp is a global operator, it manifests the 



property of essential locality []34]]. From the physical point of view, this type of models exhibit behavior 
similar to closed chains with twisted boundary conditions, however now the boundary conditions become 
sector dependent. 
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